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ABSTRACT 

We present a systematic spectral analysis with Suzaku of six AGNs detected 
in the Swift / BAT hard X-ray (15-200 keV) survey, Swift J0138. 6-4001, J0255.2- 
0011, J0350. 1-5019, J0505. 7-2348, J0601. 9-8636, and J1628. 1-5145. This is con- 
sidered to be a representative sample of new AGNs without X-ray spectral in- 
formation before the BAT survey. We find that the 0.5-200 keV spectra of these 
sources can be uniformly fit with a base model consisting of heavily absorbed 
(log Ah >23.5 cm -2 ) transmitted components, scattered lights, a reflection com- 
ponent, and an iron-K emission line. There are two distinct groups, three “new 
type” AGNs (including the two sources reported by Ueda et al. 2007) with an 
extremely small scattered fraction (/ scat < 0.5%) and strong reflection compo- 
nent ( R = 0 / 27T > 0.8 where is the solid angle of the reflector), and three 
“classical type” ones with / scat > 0.5% and R < 0.8. The spectral parameters 
suggest that the new type has an optically thick torus for Thomson scattering 
(Ah ~ 10 25 cm -2 ) with a small opening angle 6 ~ 20° viewed in a rather face-on 
geometry, while the classical type has a thin torus (Ah ~ 10 23 ~ 24 cm -2 ) with 
6 > 30°. We infer that a significant number of new type AGNs with an edge-on 
view is missing in the current all-sky hard X-ray surveys. 

Subject headings: galaxies: active — gamma rays: observations — X-rays: galaxies — 
X-rays: general 
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1. Introduction 

The growth history of supermassive black holes (SMBHs) in galactic centers is a 
key question for understanding cosmic history. This can be traced by observations of 
Active Galactic Nuclei (AGNs), where mass accretion onto SMBHs is taking place. Many 
observations indicate that obscured AGNs are the major AGN population both in the local 
universe and at intermediate-to-high redshifts (see e.g., Hasinger 2008). In fact, much of 
the SMBH growth is theoretically expected to happen during an obscured phase in the 
evolution of AGNs (e.g., Hopkins et al. 2006). Thus, detecting all the obscured AGNs in 
the universe and revealing their nature is a fundamental issue to be addressed by modern 
astrophysics. Our knowledge of this population is still quite limited even in the local 
universe, due the difficulty of detecting highly absorbed objects in most wavelength bands. 

Obscured AGNs, often referred as type 2 AGNs, are characterized by heavily absorbed 
X-ray spectra and/or by the absence of broad emission lines in the optical spectra. The 
unified model of AGNs attributes these differences to the blockage of the line of sight by a 
dusty torus. Recent X-ray and IR observations have revealed populations of AGNs without 
any signature of AGNs in the optical band (e.g., Maiolino et al. 2003). This fact implies 
that AGN surveys relying on optical emission signatures (such as the [O III] A5007 line) 
may be incomplete. 

Sensitive hard X-ray observations above 10 keV provide us with the least biased AGN 
samples in the local universe, due to the strong penetrating power against photo-electric 
absorption, except for heavily Compton thick (column density of log Ah > 24.5 cm -2 ) 
objects. Recent all sky hard X-ray surveys performed with Swift / BAT (Tucller et al. 2008) 
and INTEGRAL (Bassani et al. 2006; Krivonos et al. 2007) have started to detect many 
local AGNs with about 10 times better sensitivity than that of previous missions in this 
energy band. They have shown that roughly half of all local AGN are indeed obscured 
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AGNs. These hard X-ray samples contain a number of new AGNs, for which detailed 
follow-up studies have not been done yet. Some of these sources objects were not recognized 
to be AGNs previously even though there were optical imaging and spectroscopy data 
available. 

Our team has started to make follow-up observations of selected Swift / BAT AGNs 
with the Suzaku satellite to obtain the best quality broad band high energy data. The 
major purposes are (1) to discover new populations of AGNs, (2) to measure their spectral 
properties, which make it possible to construct the Ah distribution of a complete Swift / BAT 
sample, and (3) to accurately constrain the strength of the reflection component, which 
is an important yet uncertain parameter in population synthesis models of the X-ray 
background in explaining its broad-band shape (Ueda et al. 2003; Gilli et al. 2007). In 
particular, the simultaneous coverage over the 0.2-70 keV band with Suzaku is crucial to 
discriminate degeneracy of spectral models under possible time variability and to determine 
the amounts of scattered and reflected components, key information to constrain the 
geometry of surrounding matter around the nucleus. 

The first results obtained from the Suzaku follow-up of two Swift / BAT AGNs led to 
the discovery of a new type of deeply buried AGNs (Ueda et al. 2007). Here we present 
the results of a detailed analysis of six Swift/BAT AGNs, including the two sources (Swift 
J0138. 6-4001 and Swift J0601. 9-8636) already reported there 1 . This paper constitutes the 
first one in the series of papers to be published on Suzaku-Swift AGNs, followed by the 
second paper by Winter et al. (in preparation). § 2 summarizes the sample and observations. 
In § 3, we mainly present the results of detailed spectral analysis with Suzaku. Discussion 

1 For these two sources, the results of the present paper, which are based on more recent 
calibration of the instruments and background estimation, supersedes those by Ueda et al. 
(2007), although the essence is not changed. 
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is given in § 4. We adopt the cosmological parameters (H 0 , Q m , fl A ) = (70 km s 1 Mpc 1 , 
0.3, 0.7) throughout the paper. 


2. Observations 

We observed six Swift / BAT AGNs with Suzaku between 2006 April and 2007 January 
in the AO-1 phase. The targets are Swift J0138. 6-4001, J0255. 2-0011, J0350. 1-5019, 
J0505. 7-2348, J0601. 9-8636, and J1628. 1-5145, whose identification in the optical or 
infrared band is ESO 297-G018, NGC 1142, 2MASX J03502377-5018354, 2MASX 
J05054575-2351139, ESO 005-G004, and Mrk 1498, respectively (Tueller et al. 2008). The 
basic information of our targets is summarized in Table 1. 2MASX J03502377-5018354 and 
ESO 005-G004 had not been identified as active galaxies in the optical band, and revealed 
to contain an AGN for the Erst time with the detection of hard X-rays by the BAT. The 
sources were selected as the Suzaku targets because at the time of our proposals these 
“new” Swift AGNs had no spectral information available below 10 keV. Recently short 
follow-np observations of these targets have been carried out with X MM-Newton or the 
X-Ray Telescope (XRT) onboard Swift (Winter et al. 2008a). 

Suzaku , the fifth Japanese X-ray satellite (Mitsuda et al. 2007), carries four X-ray 
CCD cameras called the X-ray Imaging Spectrometer (XIS-0, XIS-1, XIS-2, and XIS-3; 
Koyama et al. 2007) as focal plane imager of four X-ray telescopes, and a non-imaging 
instrument called the Hard X-ray Detector (HXD; Takahashi et al. 2007) consisting of Si 
PIN photo-diodes and GSO scintillation counters. XIS-0, XIS-2, and XIS-3 are front-side 
illuminated CCDs (FI-XISs), while XIS-1 is the back-side illuminated one (BI-XIS). 

In this paper, we analyze the data of the XISs and the HXD/PIN (HXD nominal 
position), which covers the energy band of 0.2-12 keV and 10-60 keV, respectively, since 
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the fluxes of our targets above 50 keV are too faint to be detected with the HXD/GSO. 
Table 2 shows the log of the observations. The net exposure of Swift J0255. 2-0011 and 
Swift J0505. 7-2348 is about 100 ks and 80 ks, respectively, while that for the rest 4 targets 
is 20 ks each 2 . Because XIS-2 became unoperatable on 2007 November 7 (Dotani et al. 
2007), no XIS-2 data are available for Swift J0255. 2-0011 and Swift J0350. 1-5019. For the 
observations of Swift J0255. 2-0011 and Swift J0350. 1-5019, we applied spaced-row charge 
injection (SCI) for the XIS data to improve the energy resolution (Nakajima et al. 2008). 
In the spectral analysis, we also utilize the BAT spectra covering the 15-200 keV band, 
integrated for the first 9 months for Swift J0138. 6-4001 and Swift J0601. 9-8636, and for 22 
months for the rest of targets. 


3. Analysis and Results 

We analyze the Suzaku data using HEAsoft version 6.3.2 from the data products 
version 2.0 distributed by the Suzaku pipeline processing team. In extraction of the light 
curves and spectra, we set the source region as a circle around the detected position with 
a radius of 1.5-2 arcmin, depending on the flux. For the XIS data, the background was 
taken from a source-free region in the field of view with an approximately same offset angle 
from the optical axis as the source. For the HXD/PIN data, we use the so-called “tuned” 
background model provided by the HXD team. Its systematic errors are estimated to be 
~ 1.3% at a la confidence level in the 15-40 keV band for a 10 ks exposure (Mizuno et al. 


2 The primary goal of the first two targets is to constrain the reflection component with 
the best accuracy, aiming at the brightest new sources in the (then) latest BAT catalog. For 
the four targets with 20 ks exposure, the main goal was to obtain their X-ray spectra below 
10 keV to make the BAT sample complete. 
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2008). As our exposures are ~20 ks or larger, we expect that the error is even smaller than 
this value. 


3.1. Light Curves 

Figures 1 and 2 show the background-subtract light curves of our targets obtained 
with the XIS and HXD in the 2-10 keV and 15-40 keV band, respectively. To minimize 
any systematic uncertainties caused by the orbital change of the satellite, we merge data 
taken during one orbit (96 minutes) into one bin. Then, to check if there is any significant 
time variability during the observations, we perform a simple y 2 test to each light curve 
assuming a null hypothesis of a constant flux. The resultant reduced \ 2 value and the 
degree of freedom are shown in each panel. We detect no significant time variability on a 
time scale of hours for all the targets in both energy bands. Thus, we analyze the spectra 
of all the targets averaged over the whole observation epoch. 


3.2. BAT Spectra 

Before performing the spectral fit to the Suzaku data, we firstly analyze only the BAT 
spectra in the 15-200 keV band to constrain the high energy cutoff in the continuum. 

It is known that the incident photon spectra of Seyfert galaxies are well approximated 
by a power law with an exponential cutoff (cutoff power law model), represented as 
AE~ l exp{— E/E cut }, where A, T, E cut is the normalization at 1 keV, photon index, cutoff 
energy, respectively. Here we utilize the pexrav code by Magdziarz & Zdziarski (1995) to 
take into account possible contribution of the Compton reflection component from optically 
thick, cold matter. The relative intensity of the reflection component to that of the incident 
cutoff power law component is defined as R = 0/27T, where f! is the solid angle of the 



reflector (R — 1 corresponds to the reflection from a semi- infinite plane). 

In this stage, we assume R = 0 or 2 as the two extreme cases just to evaluate the 
effects of the reflection. The inclination angle is fixed at 60°. To avoid strong coupling 
between the power law slope and and cutoff energy, we fix the photon index at 1.9, which 
is a canonical value for AGNs. Table 3 summarizes the fitting results of the cutoff energy 
for each target. While it is difficult to constrain its upper limit due to the limited band 
coverage of the BAT, we find the cutoff energy must be above ^100 keV in most cases. 
Considering this limitation, we fix E cut at 300 keV in all subsequent analysis. Impacts by 
adopting a lower value of E cut in the spectral fitting are discussed in §3.5. 


3.3. Spectral Models 

We consider three basic models in the spectral analysis uniformly for all the six 
targets. Our policy is to start with the simplest model for each target; if we find that the 
fit with a simple model does not give a physically self-consistent picture and/or that the 
fit is significantly improved by introducing additional parameters, then we adopt more 
complicated models. In all the cases, we assume a power law with an exponential cutoff 
fixed at 300 keV for the incident continuum, as explained in the previous subsection. The 
Galactic absorption, Ap aI , is always included in the model (even if not explicitly mentioned 
below) by assuming the hydrogen column density from the H I map of Kalberla et al. 
(2005), available with the nh program in the HEAsoft package. For absorption, we adopt 
the photo-electric absorption cross section by Balucinska-Church & McCammon (1992) 
(“bcrnc”) and use the phabs or zphabs model in XSPEC. Solar abundances by Anders & 
Grevesse (1989) are assumed throughout the analysis. 

The first (simplest) model, designated as Model A, consists of (i) a transmitted 
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component (a cut-off power law absorbed by cold matter), (ii) a scattered component 
(a cut-off power law without absorption), and (iii) an iron-K emission line (a gaussian), 
represented as zphabs*zhighect*zpowerlw + const*zhighect*zpowerlw + zgauss 
in the XSPEC terminology 3 . According to the unified scheme, in type 2 AGNs, the incident 
power law from the nucleus will be absorbed by a dusty torus around the SMBH, while 
the nuclear emission is partially scattered into the linc-of-sight by ionized gas around the 
torus. For simplicity, we assume that the scattered component has the same shape of the 
incident power law with a fraction of / scat as first-order approximation, although, in reality, 
it often contains a number of recombination lines produced by a photo-ionized plasma (e.g., 
Sako et al. 2000) unless the scatterer is fully ionized. As a result of reprocessing by cold 
matter in the surrounding environment (such as a torus, an ionized gas, and an accretion 
disk), an iron-K emission will be produced at the rest-frame 6.4 keV. Since a line from the 
torus should not resolved by the energy resolution of the XIS, we fix the la line width of 
the iron-K line at the averaged value of the (apparent) line width of the 55 Fe calibration 
source at 5.9 keV: 0.3 eV, 38 eV, 30 eV, 0.8 eV, 1.7 eV and 3.2 eV for Swift J0138. 6-4001, 
J0255. 2-0011, J0350. 1-5019, J0505.7-2348, J0601.9-8636, and J1628. 1-5145, respectively, 
to take into account a possible systematic error in the energy response 4 . 

In the second model, Model B, we consider an additional contribution of an absorbed 
Compton reflection component from optically thick matter. This emission is expected 
from the inner wall of an (optically thick) torus and/or the accretion disk, irradiated 
by the incident continuum. Model B is expressed as zphabs*zhighect*zpowerlw + 


3 the zhighect model is utilized to represent exponential cutoff at the source redshift. 

4 The large line widths of Swift J0255. 2-0011 and J0350. 1-5019 are artifact 
caused by inaccurate calibration for the data with SCI (Matsumoto 2007; see 

http : //www . astro . isas .ac.jp/ suzaku/ analysis/xis/ver2 . 0/) 
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const*zhighect*zpowerlw + zgauss + zphabs*pexrav in the XSPEC terminology, 
where the last term represents the reflection component (not including the direct one). In 
the pexrav model, we make the solid angle hi of the reflector seen from the nucleus as a 
free parameter, and fix the inclination angle and the cutoff energy at 60° and 300 keV, 
respectively. The photon index and normalization of the power law are linked to those 
of the transmitted component. The absorption to the reflection component is set to be 
independent of that for the transmitted one, by considering a different geometry of the 
emission region. To check the physical validity of Model B, we examine the equivalent 
width (E.W.) of the iron-K emission line with respect to the reflection component, E.W. refl , 
obtained from the fit. Theoretically, it is expected to be 1-2 keV (Matt et al. 1991). This 
value, however, tightly depends on the geometries of the torus and the accretion disk. 
Thus, we regard the result of Model B as valid only if E.W. refl = 0.5 — 2 keV but discard it 
otherwise. 

We finally consider the third model (Model C), where we assume two different 
absorptions with different covering factors for the transmitted component. This reduces 
the contribution of a less absorbed reflection component compared with the case of a 
single absorption as assumed in Model B. In fact, simulations show that the absorber 
in the torus can be patchy (Wada & Norman 2002), resulting in a time-averaged 
spectrum that is better modeled by multiple absorptions to the transmitted continuum 
rather than by a single absorber. In the XSPEC terminology, Model C is expressed as 
zphabs*zpcfabs*zhighect*zpowerlw + const*zhighect*zpowerlw + zgauss + 
zphabs*pexrav, where the multiple terms of the first component (absorbed, partial 
covering model) represents the two different absorptions. 

To summarize, we can write the three models of the photon spectrum F (. E ) without 
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the Galactic absorption as follows: 

F (E) = f c exp {-N-a (E)} I (E) + (1 - / c ) exp {-A> (E)} I (E) 

+/scat/ (E) + G(E) + exp {-N^a (E)} C (E) + S (E ) , (1) 

where f c is the fraction of the more heavily absorbed component (/ c = 1 in Models A 
and B), IV-g and are the absorption column densities for the transmitted component 
(IV^ = 0 in Models A and B), a (E) is the cross section of photoelectric absorption, 

I (E) = AE~ l exp{— E/E cut } is the intrinsic cutoff power law component, / scat is the 
scattered fraction, G (E) is a gaussian representing the iron-K emission line, C (E) is the 
Compton reflection component (C (E) = 0 in Model A), and S (E) represents additional 
soft components (see below). 


3.4. Fitting Results 

Using the three models described in the previous subsection, we perform spectral fit 
simultaneously to the spectra of the FI-XISs (those of 2 or 3 XISs are summed), the BI-XIS, 
and the HXD/PIN. Based on the fitting results, we chose the best model in the following 
manner. (1) First, if Models B is found to significantly improve the fit from Model A by 
performing an F-test, we tentatively adopt Model B as a better model, otherwise Model A. 
(2) Then, if Model B is found to be physically not self-consistent in terms of the E.W. of 
the iron-K line as explained in the previous subsection, or if Model C significantly improves 
the fit compared with Model A, we adopt Model C as the best model. 

After finding the best one among the three models, we finally include the Swift /BAT 
spectrum to the Suzaku data to constrain the photon index most tightly, assuming that the 
spectral variability is negligible between the observation epoch of Suzaku and that of Swift. 
In the fit, we fix the relative normalization between the FI-XISs and the PIN at 1.1 based 
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on the calibration of Crab Nebula (Ishida et al. 2007), while those of the BI-XIS and the 
BAT against FI-XISs are treated as free parameters, considering the calibration uncertainty 
and time variability (in flux), respectively. The detailed results of spectral fit for each 
source is summarized below. 


Swift JO 138. 6-4001 

We adopt Model B as the most appropriate model of Swift J0138. 6-4001. We obtain 
(y 2 ,z/) = (101.14, 96) with Model A and (92.52, 94) with Model B, where v is the degree 
of freedom. Thus, the improvement of the fit by adding a reflection component is found to 
be significant at >90% confidence level by an F-test, which gives an F-value of 4.38 for the 
degrees of freedom of (2, 94). Since the E.W. of the iron-K line with respect to the reflection 
component is found to be E.W. rcfi > 0.5 keV, the model is physically self-consistent. Here 
we allow the absorption to the reflection component to be a free parameter, unlike the 
case in Ueda et al. (2007), where it is linked to that for the transmitted component. The 
basic results are essentially unchanged, however. 


Swift J 025 5.2-00 11 

We adopt Model B with a apec plasma model for Swift J0255. 2-0011. Models A and B 
yield (y 2 , u) = (596.85, 304) and (500.75, 302), respectively, suggesting the presence of a 
significant amount of a reflection component. We find E.W. refl = 0.81^QQg keV with Model 
B, which is physically consistent. We also confirm that applying Model C does not improve 
the fit significantly, yielding f c = 1.0 (i.e., multiple absorbers for the transmitted component 
are not required when we consider a reflection component). As positive residuals remain 
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in the energy band below 1 keV, we add the apec model 5 in XSPEC, an emission model 
from an optically-thin thermal plasma with Solar abundances, over Model B. It greatly 
improve the fit, giving (y 2 ,z/) = (330.38, 300). The temperature of the plasma is found to 
be kT = 0.741q'5q keV with an emission measure of n 2 V « 5 x 10 63 cm -3 . This temperature 
is similar to those found in many Seyfert 2 galaxies (Turner et al. 1997). 


Swift J 0350. 1-50 19 

Model C is adopted for Swift J0350. 1-5019 since we find E.W. refl < 0.1 keV . We 
obtain (y 2 , v) = (34.93, 36) with Model A, and (28.24, 35) with Model B. Here, we link 
the absorption to the reflection component and that to the transmitted one, as having 
the former as a free parameter does not help improving the fit. Model C gives (y 2 , v) = 
(28.24, 35) and E.W. refi = 2.8l];;f keV, which is now physically consistent. 


Swift JO 505.7-2348 

We adopt Model C with a line feature at ~0.9 keV as the best description for Swift 
J0505. 7-2348. We find very little reflection component in the spectra of Swift J0505. 7-2348; 
Models A and B give a very similar y 2 value of 671.3 for a degree of freedom of 602 and 600, 
respectively, with no improvement of the fit by Model B. In fact, we obtain a tight upper 
limit for the reflection component of R < 0.11 (90% confidence level) from the fit with 
Model B. Thus, we apply Model C (here A^ 00 is tied to A^ since R ~ 0). This model yields 
(y 2 , v) = (636.11, 598), which is significantly better than Model A (or B). Furthermore, 
positive residuals remain in the 0.7-1 keV band, centered at ~ 0.9 keV. To model this 


°http : // cxc . harvard . edu/atomdb/ 
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feature, we firstly add an apec component to Model C like in the case of Swift J0255. 2-0011. 
This is found to be unsuccessful, largely overestimating the 1-1.2 keV flux of the data. 
Rather, in some obscured AGNs, emission lines from Ne ions in a photo-ionized plasma 
are observed around 0.9 keV (e.g., NGC 4507 in Comastri et al. 1998 and NGC 4151 in 
Ogle et al. 2000), which is a likely origin of the excess feature seen in our spectra. As these 
emission lines cannot be resolved by CCDs, we model it with a gaussian whose la width 
is fixed at 0.05 keV, referring to the result of the Seyfert 2 galaxy NGC 4507 (Matt et al. 
2004). The center energy of the gaussian is treated as a free parameter, and is found to be 
0.87 ± 0.02keV, consistent with our assumption. We obtain (y 2 , u) = (615.08, 597), which 
is a significantly better fit compared with Model C without the line. 


Swift JO 60 1.9-863 6 

Model B is adopted for Swift J0601. 9-8636. We confirm the results reported in Ueda 
et al. (2007). Models A and B yield (y 2 ,z/) = (45.20, 37) and (32.32, 35), respectively, 
suggesting a significant amount of a reflection component. As E.W. rcfi is found to be 
2.10 ^q;3q keV, we find Model B to be a physically consistent model for this source. 


Swift J1628.1-51f5 

Model C is adopted for Swift J1628. 1-5145. Model A gives (y 2 , u) = (574.66, 546), 
which is significantly improved to (564.23, 544) in Model B but with E.W. rcfi of < 0.1 keV. 
Hence, we apply Model C to the data, which yields (y 2 , u) = (553.88, 543) and a physically 
consistent value of E.W. refl = 0.87i£| keV. 
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3.5. Results Summary 

Table 4 summarizes the choice of the spectral model (A, B, or C) and best-fit 
parameters from the Suzaku + BAT simultaneous fit for each target. The observed fluxes 
in the 2-10 keV and 10-50 keV bands and the estimated 2-10 keV intrinsic luminosity 
corrected for absorption are also listed. We find that all the six targets are heavily obscured 
with column densities larger than > 3 x 10 23 cm -2 . As already reported by Ueda et al. 
(2007), Swift J0601. 9-8636 can be called a mildly “Compton thick” AGN, since it shows log 
Ah > 24 cm -2 but a transmitted component is still seen in the hard X-ray band above 10 
keV. We confirm that the photon indices are within the range of 1. 6-2.0 for all the targets, 
and hence the fitting model is physically proper. 

The upper panels of Figure 3 show the observed spectra of the FI-XIS (black), the 
BI-XIS (red), and the HXD/PIN (magenta) folded with the detector response in units of 
counts s^ 1 keV -1 , together with the BAT spectra (blue) corrected for the detector area 
in units of photons cm -2 ks _1 keV -1 . The best-fit models are superposed by solid lines. 
The lower panels show the residuals in units of x (he., divided by the la statistical error 
in each bin). Figure 4 shows the incident spectra without Galactic absorption in units of 
EFe, where the contribution of each component in equation (1) is plotted separately; the 
black, red, blue, cyan, magenta curves correspond to the total, transmitted component, 
reflection component, scattered component, and iron-K emission line, respectively. For Swift 
J0255. 2-0011 and Swift J0505. 7-2348, the additional soft components are also included in 
purple and orange, respectively. 

For three sources Swift J0138. 6-4001, Swift J0255. 2-0011, and Swift J0601. 9-8636, 
the best fit reflection strength parameter R(= Q/2ir) > 1, although values close to unity 
are allowed within the uncertainties. Most simple geometries only allow R ~ 1 to be 
the maximum physically plausible value. The most likely explanation for these possibly 
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unphysical results is that a part of the direct emission from the nucleus is completely 
blocked by nonuniform material in the line of sight, reducing the observed normalization 
of the transmitted component smaller than the true value by a factor of < 1/R (Ueda et 
al. 2007). Alternatively, time variability can be responsible for this, if we are observing an 
echo of previously brighter phases of the AGN in the reflection component. Thus, we list 
in Table 4 the intrinsic luminosities corrected for either of these effects by multiplying R. 
Accordingly, the scattered fraction / scat is corrected by 1/R from the observed value, even 
though this correction does not affect our conclusion. 

To evaluate maximum possible errors in the spectral parameters dne to the uncertainty 
of the cutoff energy, we also perform fitting by adopting E cut = 100 keV instead of 300 keV 
for Swift J0350. 1-5019, Swift J0601. 9-8636, and Swift J1628. 1-5145, whose 90% lower limits 
of E cut are smaller than 100 keV as determined with the Swift fB AT spectra (§ 3.2). We find 
that Swift J0350. 1-5019 shows a smaller scattered fraction (/ sca t = 0.771 qu 5%) compared 
with the case of E c ut = 300 keV. For Swift J1628. 1-5145, we also obtain a slightly smaller 
scattered fraction (/ sca t = 0.74 ± 0.06) with stronger reflection ( R = 0.80^23) and larger 
absorption (IV-g = (87 ± 14) x 10 22 cm -2 ). The other parameters, including those of Swift 
J0601. 9-8636, do not change within the 90% errors. These results indicate that coupling of 
the spectral parameters with E cut could not always be negligible. Nevertheless, we confirm 
that it does not affect our discussion below, even considering the largest uncertainties as 
estimated here. 


4. Discussion 

We have for the first time obtained broad-band spectra covering the 0.5-50 keV 
band of the six “new” AGNs detected in the Swift / BAT survey that did not have precise 
X-ray observations, including the two sources already reported by Ueda et al. (2007). 
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These six targets were essentially selected without biases, and hence can be regarded as a 
representative sample of new Swift/ BAT AGNs. The spectra of all the targets are uniformly 
described with a spectral model consisting of a heavily absorbed (log Ah > 23.5 cm -2 ) 
transmitted components with a single or multiple absorptions, scattered lights, a Compton 
reflection component from optically-thick cold matter, and an iron-K emission line at 6.4 
keV in the rest-frame, with additional soft X-ray components in two cases. Thanks to its 
good sensitivity in the 10-50 keV band and simultaneous band coverage with Suzaku, we 
have accurately derived key spectral parameters including the fraction of the scattered 
component / sca t, the strength of the reflection component R(— Q/2n), and the E.W. of the 
iron-K emission line. 

The first notable result is that our sample show very small values of / scat , < 3%, 
compared with a typical value of optical selected Seyfert-2 samples of 3-10% (e.g., Turner 
et al. 1997; Cappi et al. 2006; Guainazzi et al. 2005). Winter et al. (2008b) categorize 
AGNs with / scat < 3% as “hidden” AGNs, to which all the six sources belong. This is 
not unexpected, since the “hidden” AGNs constitutes a significant fraction (24%) of the 
uniform Swift/ BAT sample. 

Figure 5 shows the correlation between the reflection strength R versus scattering 
fraction / scat for our sample. This plot suggests that there are two groups even within the 
“hidden” population defined by Winter et al. 2008b, that is, 

• R> 0.8 and / scat < 0.5% 

• R< 0.8 and / scat > 0.5%, 

although it is hard to conclude if this distribution is really distinct or continuous, given 
the small number of the current sample. Following the description by Ueda et al. (2007), 
we refer to the former group as “new type” and the other “classical type” . The new group 
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includes the two sources of Ueda et al. (2007), Swift J0138. 6-4001 and Swift J0601. 9-8636, 
and Swift J0255. 2-0011. The extremely small scattering fraction indicates two possibilities 
(1) the opening angle of the torus is small and/or (2) there is very little gas around the 
torus to scatter the X-rays. Although it is difficult to firmly exclude the latter possibility, 
the global presence of a reflection component strongly supports the former, as discussed 
below. 

We can constrain the geometry of the torus and viewing angles from the spectral 
parameters, the scattering fraction, the E.W. of the iron-K line, and absorption column 
density, based on a simple torus model. Here we utilize the calculation by Levenson et 
al. (2002), where three free parameters are introduced for a uniform-density torus whose 
cross-section has a square shape; the thickness of the torus in terms of the optical depth 
for Thomson scattering r, the half-opening angle of torus 9, and the inclination angle i. In 
this geometry, the line-of-sight column density of the torus depends on i\ when we see the 
nucleus with face-on view, the absorption becomes much smaller than the case of edge-on 
view. If the column density of the scatterer is constant, the opening angle can be connected 
to the scattering fraction as 

cos 9 = 1— j bCat (1 — cos 9 0 ) . (2) 

Jscat,0 

where we normalize / scat = / sca t,o at 9 = 9 0 . Assuming / scat ,o = 3% at 9 0 = 45° as typical 
parameters 6 , we obtain 9 ~ 20° for the new type group. Then, comparing with Figure 2 
of Levenson et al. (2002) where r = 4 is assumed, we can constrain i ~ 40° for Swift 
J0601. 9-8636, and i ~ 20° for Swift J0138. 6-4001 and Swift J0255. 2-0011, from the 
observed E.W. of the iron-K line (here we have neglected the contribution of an iron-K 
line originating from the accretion disk, which is often broader than the “narrow” one 

6 The value 9q = 45° corresponds to the number ratio of obscured AGNs to unobscured 
ones of about 3, which is consistent with the observations of local AGNs. 
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considered in the spectral model.). The absorption column densities Ah > 10 23 ' 8 cm -2 
and strong reflection R > 1 are fully consistent with the geometry considered here. The 
reflection continuum is attributable to that from the inner wall of the tall torus and that 
from the accretion disk seen with a small inclination angle. 

Thus, it is very likely that we are seeing the new type AGNs in a rather face-on 
geometry through a geometrically and optically thick torus, as discussed in Ueda et al. 
(2007). An important consequence of this interpretation is that, if we view the same objects 
in a more edge-on geometry, then these AGNs look heavily Compton-thick (Ah ~ 10 25 
cm -2 ) and hence the transmitted components hardly escape toward us even in the hard 
X-ray > 10 keV regime due to repeated scattering. Thus, simliar systems to new type AGNs 
discovered here may still be largely missed even in the on-going Swift / BAT or INTEGRAL 
surveys, requiring even more sensitive observations in hard X-rays to detect them more 
completely. 

For the classical type of AGNs, we obtain 6 > 30° from the observed scattered fraction 
by equation (2). Two sources in this group, Swift J0505. 7-2348 and Swift J1628. 1-5145, 
show weak E.W. of the iron-K emission line, ~ 0.05 keV. This is hard to explain by a 
torus with large r as assumed in Figure 2 of Levenson et al. (2002), which predicts E.W. 

> 0.1 keV if i > 6 (i.e., the transmitted component is absorbed). To consider a torus with 
different optical depths, we apply a model developed by Ghisellini et al. (1994) (see their 
Figure 1 for the definition of the geometry). We find that the E.W of the iron-K line and 
absorption can be consistently explained, within a factor of 2, if the column density of the 
torus in the equatorial plane is ~ i0 23 ~ 24 cm -2 and 60° < i < 84°. The absence of strong 
reflection components is also consistent with the small optical depth of the torus and with 
the “edge-on” accretion disk. 


To summarize, we have discovered from our first Swift- Suzaku AGN sample two 
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different classes of obscured AGN population, called “new type” that have an extremely 
small scattering fraction (/ scat < 0.5%) and a strong reflection component (R > 0.8), and 
“classical type” with / scat > 0.5% and R < 0.8. Figure 6 shows a schematic illustration of 
the torus geometry of the two types. It is likely that the new type AGNs are deeply buried 
in an optically thick torus (A% ~ 10 25 cm -2 ) with a small opening angle 6 ~ 20°, and are 
viewed in a face-on geometry. By contrast, the classical type AGNs have an optically thin 
torus (A% ~ 10 23-24 cm -2 ) with a larger opening angle 6 > 30° viewed in a more edge-on 
geometry than the new type. A significant number of new type AGNs with an edge-on 
view from us may still be missing in the current all-sky hard X-ray surveys. The presence 
of two classes of obscured AGNs implies that the understanding of AGN is not as simple 
as assumed in the population synthesis models of the cosmic X-ray background. Given the 
small number of our sample, however, it is not clear whether there are two “distinct” types 
of obscured AGNs or the distribution is continuous. More systematic follow-up observations 
with Suzaku of new Swift / BAT AGNs are important to establish the true properties of local 
AGNs. 
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^ 0 1 0 4 2x1 0 4 3x1 0 4 4x10 4 5x10 4 6x10 4 0 5x10 4 1 0 s 1.5x10 s 2x10 s 2.5x10 s 

Time (s) Time (s) 

Start Time 14062 3:40:58:184 Stop Time 14062 12:40:58:184 Start Time 13827 0:02:47:184 Stop Time 13829 3:02:47:184 



£ 0 1 0 4 2x1 0 4 3x10 4 4x10 4 5x10 4 6x10 4 o o 10 4 2x10 4 3x10 4 4x10 4 5x10 4 6x10 4 

Time(s) Time (s) 

Start Time 13838 17:58:31:184 Stop Time 13839 1:28:31:184 Start Time 13840 19:40:03:184 Stop Time 13841 9:10:03:184 


Fig. 1. — The background subtracted light curves of XIS in the 2-10 keV band during the 
Suzaku observations. One bin corresponds to 96 minutes. The data from the XIS-0 and 
XIS-2 are summed for Swift J0255. 2-0011 and Swift J0350. 1-5019, while those from XIS-0, 
XIS-2, and XIS-3 are summed for the rest. The numbers listed in each panel represent the 
value of reduced y 2 with the degree of freedom for the constant flux hypothesis. 
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(a) SWIFT J01 38.6-4001 


reduced % 2 (dot) = 1 .44 (6) 


10 4 2x1 0 4 3x10 4 4x1 0 4 5x10 4 

Time (s) 

Start Time 13890 19:02:29:184 Stop Time 13891 5:32:29:184 





(e) SWIFT J0601 .9-8636 


reduced % 2 (dot) = 0.39 (5) 


10 4 2x1 0 4 3x10 4 4x1 0 4 5x10 4 

Time (s) 

Start Time 13838 17:58:31 :184 Stop Time 13839 1 :28:31 :184 



Fig. 2. — The background subtracted light curves of HXD/PIN in the 15-40 keV band 
during the Suzaku observations. One bin corresponds to 96 minutes. The numbers listed in 
each panel represent the value of reduced y 2 with the degree of freedom for the constant flux 
hypothesis. 
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Fig. 3. — The observed spectra of the 6 targets. The black, red, magenta and blue crosses 
represent the data of the FI-XIS, BI-XIS, HXD/PIN, and BAT, respectively, showing their 
lcr errors in vertical direction. The spectra of the XIS and PIN are folded with the detector 
response in units of counts s -1 keV -1 , while those of the BAT is corrected for the detector 
area and have units of photons cm -2 ks” 1 keV -1 . The best-fit models are plotted by lines, 
and the residuals in units of % are shown in the lower panels. 
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Fig. 4. — The best fit spectral model of the 6 targets in units of EFe (where E is the energy 
and Fe is the photon spectrum). The black, dashed red, dotted blue, dot-dot-dashed cyan, 
dot-dashed magenta curves correspond to the total, transmitted one, reflection component 
, scattered component, and iron-K emission line, respectively. The purple dashed model in 
Swift J0255. 2-0011 represents the emission from an optically-thin thermal plasma, and the 
orange one in Swift J0505. 7-2348 the Ne emission lines from a photo-ionized plasma (see 
text). 




Fig. 5. — The correlation between the strength of the Compton reflection component ( R = 
0/27t) and the fraction of the scattered component (/ sca t) for the 6 targets. Filled and open 
circles represent “new type” ( R > 0.8 and / scat < 0.5%) and “classical type” ( R < 0.8 and 
/scat > 0.5%), respectively. 



Black Hole 




Classical Type 


New Type 


Classical Type 


New Type 


Fig. 6. — The schematic illustration of the torus geometry for the classical type AGNs (left 
panel) and the new type AGNs (right panel). 










Tabic 1. List of Targets 


SWIFT 

Optical/IR Identification 

R.A. (J2000) 

Dec. (J2000) 

Redshift 

Classification 

J0138. 6-4001 

ESO 297-G018 

01 38 37.16 

-40 00 41.1 

0.0252 

Seyfert 2 

J0255. 2-0011 

NGC 1142 

02 55 12.196 

-00 11 0.81 

0.0288 

Seyfert 2 

J0350. 1-5019 

2MASX J03502377-5018354 

03 50 23.77 

-50 18 35.7 

0.036 

Galaxy 

J0505. 7-2348 

2MASX J05054575-2351139 

05 05 45.73 

-23 51 14.0 

0.0350 

Seyfert 2 

J0601. 9-8636 

ESO 005-G004 

06 05 41.63 

-86 37 54.7 

0.0062 

Galaxy 

J1628. 1+5145 

Mrk 1498 

16 28 4.06 

+51 46 31.4 

0.0547 

Seyfert 1.9 


Note. — The position, redshift, and classification for each source is taken from the NASA/IPAC Extragalactic 
Database. 


Table 2. Observation Log 


SWIFT 

Start Time (UT) 

End Time 

Exposure 3, (XIS) 

Exposure (HXD/PIN) 

SCI b 

J0138. 6-4001 

2006 Jun 04 18:13 

Jun 05 05:00 

21.2 ks 

16.5 ks 

Off 

J0255. 2-0011 

2007 Jan 23 14:54 

Jan 26 05:30 

101.6 ks 

80.6 ks 

On 

J0350. 1-5019 

2006 Nov 23 02:07 

Nov 23 13:33 

19.3 ks 

15.7 ks 

On 

J0505. 7-2348 

2006 Apr 01 22:12 

Apr 04 02:55 

78.6 ks 

50.3 ks 

Off 

J0601. 9-8636 

2006 Apr 13 16:24 

Apr 14 01:52 

19.8 ks 

15.7 ks 

Off 

J1628. 1+5145 

2006 Apr 15 18:49 

Apr 16 09:22 

23.6 ks 

19.5 ks 

Off 


a Based on the good time interval for XIS-0. 

b With/without the spaced-row charge injection for the XIS(Nakajima et al. 2008). 
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Tablc 3. 


Cutoff Energies (E cut ) determined by the BAT spectra 


Cl/ 2n 

SWIFT J0138. 6-4001 

J0255. 2-0011 

J0350. 1-5019 

J0505. 7-2348 

J0601. 9-8636 

J1628. 1+5145 

0 

> 317 

> 289 

> 106 

> 318 

> 116 

> 99 

X 2 /d.o.f. 

1.58/2 

4.57/6 

1.16/6 

13.58/6 

4.96/2 

2.38/6 

2 

329 (> 91) 

> 253 

> 95 

> 306 

> 82 

> 87 

X 2 /d.o.f. 

1.05/2 

5.21/6 

2.40/6 

15.24/6 

3.37/2 

4.61/6 


Note. 


The unit of _E C ut is keV. 



Tabic 4. Best-fit Spectral Parameters 
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transmitted one in Model C. (5) The power-law photon index. (6) The fraction of the scattered component relative to the intrinsic power law, corrected 
for the transmission efficiency of 1/R if R > 1 (see text). (7) The center energy of the iron-K emission line at the rest frame of the source redshift. (8) 
The observed equivalent width of the iron-K line with respect to the whole continuum. (9) The line-of-sight hydrogen column density for the reflection 
component. (10) The relative strength of the reflection component to the transmitted one, defined as R = S7/27T, where 17 is the solid angle of the 
reflector viewed from the nucleus. (11) The observed flux in the 2-10 keV band. (12) The observed flux in the 10-50 keV band. (13) The 2-10 keV 
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